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We investigate the bounds on axionlike states from flavor-changing neutral current b — > s decays, 
assuming the axion couples to the standard model through mixing with the Higgs sector. Such 
GeV-scale axions have received renewed attention in connection with observed cosmic ray excesses. 
We find that existing B — > K£ + £~ data impose stringent bounds on the axion decay constant in the 
multi-TeV range, relevant for constraining the "axion portal" model of dark matter. Such bounds 
also constrain light Higgs scenarios in the next-to-minimal supersymmetric standard model. These 
bounds can be improved by dedicated searches in B-factory data and at LHCb. 
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I. INTRODUCTION 

Motivated by a variety of cosmic ray anomalies , 
a new dark matter paradigm has emerged where TeV- 
scale dark matter interacts with GeV-scale bosons 
In one such scenario — dubbed the "axion portal" 
- dark matter in the Milky Way halo annihilates into 
light pseudoscalar "axions" . In order to explain the ob- 
served galactic electron/positron excess, the axion, a, is 
predicted to have a specific mass and decay constant B 

360 < m a < 800 MeV, f a ~ 1 - 3 TeV. (1) 

These axions couple to standard model fermions propor- 
tional to their Yukawa couplings, and in this mass range 
the axion dominantly decays as a —> /i + fi~ . Other novel 
dark matter scenarios involving axionlike states have also 
been proposed p|-|l3|], which allow for broader range of 
axion masses and decay constants. 

More generally, light axionlike states appear in a va- 
riety of new physics scenarios, as they are the ubiqui- 
tous prediction of spontaneous Peccei-Quinn (PQ) Q 
symmetry breaking. The most famous example is the 
Weinberg- Wilczek axion invoked to solve the strong CP 
problem ]l5| , as well as invisible axion variants |l7| - 
|20f . Light pseudoscalar particles appear in any Higgs 
sector with an approximate PQ symmetry, which often 
occurs in the minimal or next-to-minimal supersymmet- 
ric standard models (MSSM and NMSSM). Models of 
dynamical supersymmetry breaking typically predict an 
_R-axion [ pj| , whose couplings can mimic PQ-type axions. 
There has also been speculation [^2| that the HyperCP 
anomaly p3| might be explained by a light axion. There- 
fore, searches for light axionlike states have the potential 
to confirm or exclude a variety of new physics models. 

In this paper, we show that flavor-changing neutral 
current b —> s decays place stringent bounds on such 
models. While the coupling of the axion to fermions is 
flavor-diagonal, the b —> sa decay mediated by a top-IF 
penguin diagram is enhanced by the top Yukawa coupling 
appearing in the top-axion vertex. To our knowledge, 
Refs. [p4| p6[ were the first to consider this decay as a 
search channel for light pseudoscalars, where the a field 
was identified with the CP-odd Higgs A in a two Higgs 
doublet model (2HDM). The goal of this paper is to re- 



vive this search channel in models like the axion portal, 
where there is an a field which mixes with A . 

In the parameter range of interest for the axion portal, 
the axion decays promptly to and we find that 

existing B — > K£ + £~ data (for £ — e, \i) can be used to 
derive multi-TeV constraints on the axion decay constant 
f a , especially for small values of tan/3. For heavier ax- 
ionlike states with reduced branching fractions to muons, 
B K£ + £~ can still be used to place a bound, relevant 
for constraining light Higgs scenarios in the NMSSM If- 
p^ |. The estimates in this paper are likely improvable 
by dedicated B — > Ka searches at BaBar and Belle, 
and can be further strengthened at LHCb and a possible 
super £?-factory. These searches are complementary to 
T(nS) 7a searches recently performed by BaBar [ |30| . 

In the next section, we review the axion portal La- 
grangian, which is relevant for any DFSZ-type (Dine- 
Fischler-Srednicki-Zhitnitsky) axion jl9[ po| , and use it 
to calculate the effective b — > sa vertex in Sec. Ej. We 
sketch the current experimental situation in Sec. JV| and 
derive corresponding bounds in Sec. [v[ We conclude in 
Sec. [V3. 



II. REVIEW OF THE AXION PORTAL 

If one were only interested in studying the tree-level 
interactions of new axionlike states, it would be sufficient 
to introduce a new term in the Lagrangian of the form 



Ja 



(2) 



where f a is the axion decay constant and is the fermion 
charge under the broken 17(1). By the equations of 
motion, such a coupling is proportional to the fermion 
mass parameter, leading to an effective coupling constant 
tym^ I f a . However, the b — > sa process we are interested 
in occurs via a top- IF penguin loop. With only Eq. (Q), 
such a diagram is logarithmically sensitive to the cutoff 
scale [53] , so it is necessary to embed the axion coupling 
in a complete theory to get a reliable bound on f a . 

The axion portal || is an example of a class of theories 
where the b — > sa amplitude is finite. The axion arises 
from spontaneous PQ-symmetry breaking in a 2HDM, of 
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which the DFSZ axion is a special case. We show that 
the b — > sa amplitude can be derived from the b — > sA° 
amplitude, where A is the CP-odd Higgs boson in a 
PQ-symmetric 2HDM. 

Consider a complex scalar field S carrying J7(1)pq 
charge that gets a vacuum expectation value (S) = f a - 
This spontaneous symmetry breaking leads to a light ax- 
ionlike state, a, 



S = f a exp 



V2f a 



(3) 



The assumption in the axion portal (and for any DFSZ- 
type axion) is that the only operator that transmits PQ 
charge from S to the standard model is 



SC = \S n h u h d + h.c. 



(4) 



where A is a (possibly dimensionful) parameter, and n is 
an integer. This coupling forces h u hd to carry nontrivial 
PQ charge, and we assume that the entire Higgs potential 
preserves the £/(1)pq symmetry to a good approxima- 
te DFSZ axion || " 
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corresponds to the case 
while for the PQ-symmetric NMSSM pi 



tion. 

with n = 2, 

n = 1. Either case can be used in the axion portal model 
of dark matter. 

Since the vacuum expectation values (VEVs) of S, h u , 
and hd all break the PQ symmetry, the physical axion will 
be a linear combination of the phases of all three fields. 1 
However, in the /„ ^> uew limit, it is calculationally 
more convenient to work in an "interaction eigenstate" 
basis, where the axion a only appears in S, and the CP- 
odd Higgs A only appears in the two Higgs doublets in 
the form: 



h u = 



hd 



v u exp 



i cot f3 a0 







Vd exp 





i tan j3 
%/2-uew 



where tan j3 = v u /vd, 



vew 



9 



174 GeV, 



(5) 



(6) 



and we have omitted the charged Higgs P ± and the CP- 
even Higgses for simplicity. The coefficients appearing in 
front of ^4° ensure that ^4° is orthogonal to the Goldstone 
boson eaten by the Z boson. 

This exponential parameterization of ^4° is convenient 
for our purposes, since PQ symmetry implies that mass 



terms involving a and A can only appear in Eq. (0). In 
this basis, the physical degrees of freedom are given by 



Ophys. = a cos 9 — A sin 9, 
A° hys = a sin 8 + A cos 9, 



with 



tan# 



uew sin 2/3 
17 2 



(7) 



(8) 



At this level, the physical axion is massless. 2 A small 
mass (beyond the contribution from the QCD anomaly) 
can be generated by a small explicit violation of the PQ 
symmetry, but the precise way this happens is irrelevant 
for our discussion. 

The dominant decay mode for the axion depends on 
its mass, m a . The axion decay width to an £ + £~ lepton 
pair is given by 



16tt fi 



4m? 



(9) 



For 2m e < m a < Im^, the dominant decay is a — > e + e~ . 
In this mass range, however, strong bounds already exist 
from K — > ira decays |52[ |33|. With the axion decay to 
fermions being proportional to their mass-squared, a — > 
dominates over a — > e + e~ for m a > 2m n . Note 
that in the mass range given in Eq. ([!]), the axion decays 
within the detector as long as /„ < 1000 TeV. 

The axion decay becomes more complicated at higher 
masses when hadronic decay modes open up. Refer- 
ence H estimated that the a — > 3ir channel starts to 
dominate over the channel at m a ~ 800 MeV. 

Hadronic channels dominate the axion decay until m a > 
2m T , when the t + t~ channel becomes dominant. How- 
ever, as emphasized recently in |p4|| , throughout the en- 
tire mass range 2m n < m a < 2m;,, the branching ratio 
to remains significant, and until the t + t~ thresh- 

old, it never drops below 0(1O -2 ). For m a > 2m T , the 
branching fraction to muons is approximately 



Br(a — > 



0.003, 



(10) 



with the precise value depending on tan f3 through T(a 
cc) and on the neglected phase space factor. 



III. THE EFFECTIVE b -> sa COUPLING 

By assumption, the physical axion state dominantly 
couples to standard model fields via its mixing with A . 



1 This also means that the physical axion decay constant will be 
a function of the three VEVs. The difference is negligible when 
fa S> dew i and we will continue to refer to f a as the axion decay 
constant. 



2 For completeness, the physical A mass is given by m 2 (A® h y s ) = 
A(/ a )"(2/sin2/3)(l + tan 2 6»). 
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Therefore, at one-loop level, the amplitude for b sa 
can be derived from 



Ml 



sm9 x M(b -»■ sA°) 2IiUM , (11) 



where "2HDM" refers to a (PQ-symmetric) 2HDM with 
no S field. Moreover, since the final state only contains 
a single axion field, there is no difference in the relevant 
Feynman rules between the exponential parameterization 
in Eq. @ and the standard linear parameterization of A° 
in the two Higgs doublet literature. For concreteness, we 
will consider a type-II (MSSM-like) 2HDM. 3 

The radiatively induced b — > sA° coupling in a type-II 
2HDM was calculated in the early 1980's independently 

M 



m two papers 



2o 



The dominant contributions come 
from penguin diagrams involving a top quark, a W bo- 
son and/or charged Higgs boson, and the ttA° or 
W^H^A couplings (and corresponding counterterms). 
The one-loop b — > sA° amplitude is reproduced to lowest 
order (in the mB,A° *C mw,t,H limit) by the tree-level 
matrix element of the effective Hamiltonian [B9, Bql 4 



n = 



fV t * s V tb m\ 



128 7T 2 



(X x cot p + X 2 cot 3 P) s^P L b df,A°. 

(12) 

The functions Xi and X 2 depend on the charged Higgs 
boson mass to#, and are given by 



Xi = 2 



3m^, 



m H — — m 



w 



3mVK + m w- 2m ?) ln m t 



[m 2 H - rrtfy) (m 2 - m^) 2 m 2 w 



'H 



H 



'H 



H 



mi m 



H 



w 



ln-o-, 



•H 



x 9 = - 



2m 2 



9 2 

m H — m\ 



2 2 

mi, , mi 
- In — - 

n 2 ™2 ™2 



m 



H 



(13) 



From this effective Hamiltonian, we can calculate vari- 
ous B decay rates in the 2HDM. These are summarized in 
the Appendix for B — > Ka, B — > K*a, and the inclusive 
B X s a rates. Using Eq. (flT|), the rates in any of these 
channels relevant for the axion portal arc determined by 



T(B Ka) = sin 2 x T(B 



KA°) 2R u M . (14) 



The type-I 2HDM model gives the sam p b — > sA° amplitude to 
the order we are working; see Ref . fca] . 

The results published in these two papers differ, a fact which 
seems to have gone unnoticed — or at least unremarked upon — 
in the literature. We have redone the calculations both in the 
unitary gauge and in the Feynman gauge and agree with the 
result in Ref. pq | . We also agree with Ref. p6| if we replace in 
their Eq. (9) the second ln(r7ij/m^) term by \n{mf / m 2 H ) , most 
likely indicating a simple typographical error. Severa l pa pers in 
the literature seem to use the result as printed in Ref. ]26| , which 
has qualitatively wrong implications. For example, it exhibits 
decoupling in the ran ~> 00 limit and singularities when ran ~* 
mt, whereas the correct result does not. 



IV. EXPERIMENTAL BOUNDS 

In the parameter range of interest, the axion has a 
significant decay rate to leptons and decays promptly on 
collider timescales. Thus, the axion would manifest itself 
as a narrow dilepton peak in b — > s£ + £~ decays. 

The b — > sa — > s £ + £~ process contributes to both in- 
clusive and exclusive B — > X s £ + £~ decays [j3l| |36f . These 
final states receive large long-distance contributions from 
intermediate J/ip and ip' resonances decaying to £ + £~, 



which result in removing the surrounding q {= 



■) 



regions from the measurements. The so-called low-g 2 re- 
gion (q 2 < 7 — 8 GeV 2 ) can probe axion masses up to 
m a ~ 2.7 GeV, while the high-g 2 region (q 2 > 14 GeV 2 ) 
is above the a — > t + t~ threshold. In general, one can 
bound the axion contribution in both these regions. 

In the low-g 2 region, and especially for m a < 800 MeV 
as in Eq. (|l|), the exclusive mode B — > K£ + £~ is par- 
ticularly well-suited to constrain b — > sa. This is be- 
cause dT(B — > K£ + £~)/dq 2 varies slowly at small q 2 , and 
B -S- K£+£" has a smaller rate than B -> K*£+£-, thus 
it gives us the best bound by simply looking at the mea- 
sured spectrum. In contrast, the exclusive B —> K*£ + £~ 
and the inclusive B — > X s £ + £~ decay modes receive large 
enhancements from the electromagnetic penguin opera- 
tor, whose contribution rises steeply at small q 2 , as 1/q 2 . 
This will complicate looking for a small excess in these 
modes in this region. 

For m a > 1 GeV, we expect that the bounds from 
B — > K£ + £~ and K*£ + £~ may be comparable (possi- 
bly even from B — > X s £ + £~ if a super B-factory is con- 
structed), and a dedicated experimental analysis should 
explore how to set the strongest bound, using the rate 
predictions in App. . For the remainder of this paper, we 
focus on B -» K£+£-. 

Since B —> Ka contributes mostly to the K/i + yT final 
state, and much less to Ke + e~, to set the best possible 
bound on B — > Ka, one needs the B — > Ku, + ur and 
B — > Ke + e~ spectra separately. This information does 
not seem to be available in the published papers p?], |8| . 
Based on the latest world average, Br(£> — > Kv~£~) — 
(4.5 ± 0.4) x 10~ 7 (37|-(39|, and the spectrum in Fig. 1 in 
Ref. Wi\, it seems to us that 



Br(S Ka) x Br(a ^ + /0 < 10" 



(15) 



is a conservative upper bound for any value of the axion 
mass satisfying m a < mB — mx- 

As we emphasized, BaBar, Belle, and a possible super 
B-factory should be able to set a better bound on a nar- 
row resonance contributing to B — > K^*> /i + u~ but not 
to B -)• K^e+e-. Moreover, LHCb will also be able 
to search for deviations from the standard model predic- 
tions in B — > K(*H + £~ with significantly improved sen- 
sitivity. While we could not find a recent LHCb study 
for the K£+£- mode (only for K* £+£-), the fact that 
the signal to background ratio at the e + e~ 23-factories is 
not worse in B K£ + £~ than in B ->• K*£ + £r suggests 
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Bound on f a tan /? (Large tan /?) 



FIG. 1: Bounds on f a as a function of tan /3 and mn for n — 1 
in Eq. (^), for m a <C rrij. For each displayed value of f a there 
are two contour lines, and the region between them is allowed 
for f a below the shown value. The bound disappears along 
the dashed curve, and gets generically weaker for larger tan /3. 



that LHCb should be able to carry out a precise mea- 
surement jh|. Interestingly, since the B — > Ka signal is 
essentially a delta function in q 2 , the bound in Eq. ( flo] ) 
can be improved as experimental statistics increase by 
considering smaller and smaller bin sizes, without being 
limited by theoretical uncertainties in form factors |4l|| 
(or by nonperturbative contributions The bound 

on f a will increase compared to the results we obtain in 
the next section, simply by scaling with the bound on 
l/V / BrXB z > Ka). 



V. INTERPRETATION 

We now derive the bounds on f a using the calculated 
B — > Ka branching ratio in Eq. (|l4j) and the experimen- 
tal bound in Eq. (|l5|). We start with the axion portal 
scenario with Br(et — > /i + /i~) ~ 100% and where sinfl is 
defined in terms of f a by Eq. (||). We will then look at 
the bound on more general scenarios, including the light 
Higgs scenario in the NMSSM. 

For the axion portal, Fig. [j] shows the constraints on f a 
as a function of the charged Higgs boson mass mn and 
tan/3. For concreteness, we take n = 1; other values of n 
correspond to a trivial scaling of f a . In the mass range 
in Eq. ([j]) , the dependence on m a is negligible for setting 
a bound. The bound on f a is in the multi-TeV range for 
low values of tan/3 and weakens as tan/3 increases. At 
each value of tan /3, there is a value of mu for which the 
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large tan/3 limit. To indicate the region of validity of the 
large tan/3 approximation, the dashed (dotted) curve shows 
the bound for tan/3 = 3 (tan/3 = 1). 



b — > sa amplitude in Eq. (12[) changes signs, indicated by 
the dashed curve in Fig. l|7 along which the bound dis- 
appears. Higher order corrections will affect where this 
cancellation takes place, but away from a very narrow re- 
gion near this dashed curve, the derived bound is robust. 
The region tan /3 < 1 is constrained by the top Yukawa 
coupling becoming increasingly nonpertubative; this re- 
gion is included in Figs. [l] and nevertheless, to provide 
a clearer illustration of the parametric dependence of the 
bounds. 

As one goes to large values of tan /3, the X\ piece 
of Eq. (fij) dominates, and sin(2^)/2 = l/tan/3 + 
0(1/ tan /!3). In this limit, the constraint takes a par- 
ticularly simple form that only depends on the combi- 
nation f a tan 2 P, as shown in Fig. |^. Except in the re- 
gion close to tuh ~ 550 GeV, the bound is better than 
f a tan 2 /3 > few x 10 TeV. 

These B — > Ka bounds are complementary to those 
recently set by BaBar J3(| in Y(nS) — >• ja — >• 7/i + /i _ : 



f a > (1.4 TeV) x sin 2 /3. 



(16) 



For example, for mn — 400 GeV, the T bound dominates 
for tan/3 > 5, while B — > Ka dominates for tan /3 < 5. 

The bounds in Figs. |l| and ^| apply for a generic axion 
portal model where mu and tan /3 are free parameters. 
One would like some sense of what the expected values 
of mn and tan/3 might be in a realistic model. Ref. || 
considered a specific scenario based on the PQ-symmetric 
NMSSM Hi). In that model small tan ,3 is preferred, 



since large tan /3 requires fine-tuning the Higgs potential. 
In addition, mn is no longer a free parameter and is 
approximately related to the mass of the lightest CP- 
even scalar sq via 



2 m So fg 
sin 2 2/3 dew 



(17) 
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FIG. 3: Bounds on sin 2 #Br(a as a function of 

tan/3 and tuh- Similar to Fig. [I], the successively darker 
regions going away from the dashed curve are allowed for 
sin 2 9 Br(a — > fi + fi~) above the indicated values. When m a is 
not small compared to mg, these bounds should be modified 
by Eq. (H), but this is a small effect. 

In the context of dark matter, Ref. [§J required m SQ to be 
O(10 GeV) to achieve a Sommerfeld enhancement. Tak- 
ing m So = 10 GeV and f a = 2 TeV as a benchmark, the 
B — » Ka bound requires 2.5 < tan/3 < 3.0, correspond- 
ing to 490 GeV < m H < 650 GeV, in the vicinity of 
the cancellation region. This bound is very sensitive to 
to So ; for m So = 20 GeV and f a = 2 TeV, the bounds are 
1.5 < tan ,3 < 1.7 and 550 GeV < m H < 610 GeV. Note 
that models like JnJ have no preferred value of m# , can 
have larger values of f ai and do not disfavor large tan /3. 

As mentioned, these B — ¥ Ka constraints apply to any 
scenario where the branching ratio formula in Eq. ( |T4| ) 
applies, i.e. where the axion couplings are determined via 
Eq. ^ , and where m a < raj — wlk ■ For example, recent 
studies of light Higgs bosons in the NMSSM |27|-|f] and 
related dark matter constructions [[l0| |l2| also contain 
a light pseudoscalar whose couplings to standard model 
fermions can be described in terms of a mixing angle 9, 
as m Eq. (§). 5 There, the mass of the a field is expected 
to be 2m T < m a < 2m{,, with the a — > /i + /i~ branching 
fraction estimated in Eq. (|lO|) . 

To show the constraints on such scenarios in a model 
independent way, we plot the bound on the combination 



5 In the literature, sinfl is often referred to as the "non-singlet 
fraction" cos 8 a ¥■% ■ 




m H (GeV) 

FIG. 4: Bounds on sin 2 #Br(a — } fi + fj,~)/ tan 2 (3 in the large 
tan/3 limit. The shaded region is excluded, and the dashed 
(dotted) curve shows tan /3 = 3 (tan/3 = 1). 

sin 2 9 Br(a — > in Fig. |], in the m 2 , <?C m 2 B limit for 

simplicity. We also show the large tan /? limit in Fig. |^, 
where the bound is on the combination sin 2 9 Br(a — >■ 
tan 2 /?. To apply these bounds for the case where 
m a is not small compared to mj, one should make the 
replacement in Figs. || and [| (see the Appendix), 

■2/}^ • 2 q X K(m a )[fo(m 2 a )] 2 . 2 
sin 9 sin 9 ^esih 9 R(m a ) . 

(18) 

Using a simple pole form for the q 2 dependence of fo [pTJ , 
we find that R(m a ) deviates from unity by less than 20% 
for m a < 4.6 GeV (i.e. nearly over the full kinemati- 
cally allowed region), and so it is a good approximation 
to neglect R(m a ). In the case of NMSSM scenarios, the 
precise bound depends strongly on the parameters of the 
theory. To give a sense of the strength of the bound, for 
m a ~ 4 GeV, m H - 200 GeV, and using Eq. (|l0|), the 
bound at large tan/3 implies sin 2 9/ tan 2 < 5 x 10~ 4 
(and sin 2 9 < 2 x 10~ 4 for tan /3 = 1), which is a signifi- 
cant constraint on large mixing angles or small tan f3. 



VI. CONCLUSIONS 

In this paper, we explored bounds on axionlike states 
from flavor-changing neutral current b — > s decays. We 
found that the exclusive B — > K£ + £~ decay is particu- 
larly well-suited to constrain such contributions. In the 
case of the axion portal (or equivalently, any DFSZ-type 
axion), we derived a bound from current -B-factory data 
on the axion decay constant f a . The bound is in the 
multi-TeV range, gets stronger for small tan/3, and de- 
pends sensitively on the value of the charged Higgs bo- 
son mass. This places tension on the axion portal model 
of dark matter in the parameter space given in Eq. (Q). 
More generally, there is a constraint on any pseudoscalar 
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with 2m ^ < m a < tub — tuk whose couplings to standard 
model fermions arise via mixing with the CP-odd Higgs 
A . This is true even if Br(a -> n + [J.~) ~ O(10~ 3 ), as is 
the case for light Higgs scenarios in the NMSSM. 

We derived our bound using a conservative estimate 
from the q 2 distribution in B — > K£ + £~ . The bound 
could most probably be improved through a dedicated 
search in existing B-factory data, and in searches at 
LHCb and a possible future super £>-factory. The 
B — > Ka search is complementary to axion searches in 
T(nS) — > 7a, because for fixed mixing angle 9 in a type- 
II 2HDM, the former scales like 1/ tan 2 f3 while the latter 
scales like tan 2 (3. 

One way to extend our analysis would be to look at 
axions decaying to hadronic final states. We focused on 



the decay mode a 



since the a — > e + e mode 



is already well-constrained by kaon decays, and we were 
motivated by the parameter space relevant for Ref. [||. 
However, as the axion mass increases, other decay chan- 
nels open up, such as a — > n + n~TT°, a — > KK*, etc. 
These would also be worthwhile to search for in £?-factory 
data, especially since dark matter models such as IllJ arc 
compatible with a — > 7r + 7r~7r° decays. It appears to us 
that setting bounds in these modes is more complicated 
than for B — > K£ + £~, and should be done in dedicated 
experimental analyses. For constraining higher mass ax- 
ions, it would be interesting to study whether B-factories 
could search for narrow resonances in B — > Kt + t~ at a 
level of sensitivity no weaker than m 2 / m 2 times the cor- 
responding bound in B — > Kfi fi~. Combining a num- 
ber of search channels, one would be able to substantially 
probe scenarios containing light axion-like states. 

Note Added: While this paper was being completed, 
Ref. ifill appeared, which claims much stronger bounds 
on f a than our result. They use a different effective 
Hamiltonian from Eq. ([jj]) , which does not include the ef- 
fect of charged Higgs bosons, crucial for bounding DFSZ- 
type axions. 
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Appendix: Decay Rates 

In this Appendix, we list the B decay rates to K^*>A° 
and X S A° in the 2HDM, using the effective Hamiltonian 



in Eq. . These should be combined with Eq. (|l4|) to 
bound the axion models. 
Defining 

r = G ^J^ tb \ 2 m* m \{X x cot /3 + X 2 cot 3 p) 2 , (A.l) 



V2 2 12 tt 5 



and 



X K( ,) = \J (m| - m 2 A0 - m 2 KM ) 2 - 4m 2 A0 m 2 KM 



the B —> KA° decay rate is given by 



T(B ->■ KA a ) = T a 



\ K {m 2 B -m 2 K ) 2 



(A.2) 



[/oK»f (A.3) 



The B -> K*a decay rate is 



r(B^r/) = r %[4 (™^)] 2 . (A.4) 



In both decays we used the standard definitions |3E| of 
the form factors, 

(K(p-q)\s4P L b\B(p)) - \{m 2 B -m 2 K )f Q {q 2 ), (A.5) 
(K*(p - q)\HP L b \B(p)) - -i m K , (e* ■ p) A (q 2 ) . 

(We caution the reader not to confuse ^4° and A , each 
of which are standard in the respective contexts.) 

In Eq. (A.l) it is the MS top quark mass which enters, 
appropriate both for the coupling to Higgses and in loop 
integrals. While this distinction is formally a higher order 
correction, since the rates are proportional to mf , we use 
the Tevatron average top mass, converted to MS at one- 
loop, m t — m t [l — 4a s /(37r)] R3 165 GeV. 

The largest hadronic uncertainty in evaluating the im- 
plication of the bound in Eq. ( p"5| ) is the model depen- 
dence in the calculations of the form factor / (to 2 ), which 
is an increasing function of q 2 . For /o(0), QCD sum 
rule calculations obtain values around 0.33, with an or- 
der 10% uncertainty [Q. To be conservative, in eval- 
uating the bound on / a , we only assume /o(0) > 0.25 
for m a <^ tub (which also covers lower values moti- 
vated by the fit in Ref. Q]). For m a > 2m r , rele- 
vant for Eq. (113), we use the approximation fo(q 2 ) = 



/ (0)/(1 - g 2 /37.5 GeV 2 ) which should be good 
enough for our purposes. For recent QCD sum rule cal- 
culations of Ao(q 2 ), relevant for setting a bound using 
B -» K*£+£-, see Ref. g§. 

The inclusive B — > X s a decay rate, which can be cal- 
culated (strong interaction) model independently in an 
operator product expansion, is given at leading order in 
^QCu/m b by 



T(B^X S A°)=2T ^I- 1 



A 



(A.6) 
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